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^_2 Abstract. Theta graphs are important geometric graphs that have many applications^ 

(^ including wireless networking, motion planning, real-time animation, and minimum- 

CN spanning tree construction. We give closed form expressions for the average degree of 

theta graphs of a homogeneous Poisson point process over the plane. We then show that 
essentially the same bounds — with vanishing error terms — hold for theta graphs of finite 
sets of points that are uniformly distributed in a square. Finally, we show that the num- 
ber of edges in a theta graph of points uniformly distributed in a square is concentrated 
around its expected value. 
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Figure 1: A vertex, u, in a Sy-graph; (a) u connects to the "nearest" point in each cone; 
(b) "nearest" is defined in terms of the projection onto the axis of each cone. 



1 Introduction 

Theta graphs [1101 [T6l 117] are important geometric graphs that have many applications, 
including wireless networking ^, motion planning [IIOIJ , real-time animation [14], and 
minimum-spanning tree construction [I24| |. These graphs are defined on a planar point set, 
S, with an integer parameter k. For each i e{l,...,k}, define the cone 

Q = {u 6 R^ : Iqou 6 [2n{i - l)/k,2ni/k)} , 

where q = (1,0) and - (0,0). In the 0;t"graph, 0yt(S), each point u ^S has an edge connect- 
ing it to the nearest point, if any, in the cone Ci+u, for each i € {!,..., k). Here, "nearest" has 
a special meaning: The theta graph connects u to the point whose orthogonal projection 
on the axis of C, -h u is closest to u (see Figure [TJ. 

Theta graphs have two important properties that make them suited to a wide va- 
riety of applications: They are sparse; d]^{S) has at most k\S\ edges and they are spanners; 
the length of the shortest path between any two vertices u and w is at most a constant 
(depending only on k and not on S) times the Euclidean distance between u and w. For 
any point set, S, di{S) is a spanner for any fc > 4 [f5l[7l[8l [T6ll23] . 

Note that, although 0jt(S) has at most k\S\ edges, it can also have significantly fewer 
edges. For example, if the points of S all lie on a line, then 0jt(S) has only \S\ - 1 edges. 
More typical, though, is for 0yt(S) to have somewhere between k\S\/2 and k\S\ edges; each 
vertex u eS chooses k edge^of the graph but sometimes an edge uw is chosen both by u 
and w and therefore it should not be counted twice. 



Specifically only the points that are unoriented 6-maxima of S may choose fewer than k edges; the number 
of unoriented 0-maxima is typically O(VfSl) I4j Theorem 4]. 



1.1 The Models and Results 

In this paper, we study the typical number of edges in dj^-graphs by studying the average 
number of edges in two different models of random point sets. 

We begin, in Section [2j by studying (the infinite) 0jt -graphs generated by a homo- 
geneous Poisson point process with unit intensity over the entire Euclidean plane. In this 
model, we study the average degree of a vertex in the 0/;-graph. For the 0jt"graph^ this 
quantity is at most 2k since each vertex defines k edges of the graph and each edge has 
2 endpoints. However, in some cases an edge uw is mutual in the sense that the edge is 
created both by u and by w. If we let p/f denote the probability that an edge of the 0;t"gr^ph 
is mutual, then the average degree of the 0j. graph is 

(The second term corrects for the double-counting of mutual edges.) Thus, understanding 
the average degree of a Oj^-graph. boils down to computing pi^. 

In Section [2?T] we show that, for all even integers k> A, 

7tV3 
Pk^ -^ « 0.6045997883 . 

That is, the probability of an edge being mutual is independent of k. Thus, for all even 
integers fc > 4, the average degree of a vertex is 

dk^\2-——]kxl.395A002l2-k . 



In Section [2?2] we show that, for odd values of A: > 5, the situation is very different. 
The mutual edge probability, p^, depends on fc in a complicated way that includes trigono- 
metric functions and square roots. However, the value of p^ is significantly larger than 

^^ for all odd values of k. Indeed, p^ is a decreasing function of k and 
p^ > lim pk = 2arctan(l/3) « 0.6435011088 . 

fc— >oo 

Thus, for all odd values of fc > 5, 

dk<(2- 2arctan(l/3))fc « 1.356498891/: 
Thus, in some sense, odd values of k offer "more bang for the buck." 

In Section[3l we also study the i.u.d. model, in which a set, S, of n points is indepen- 
dently and uniformly distributed in a square. In this model, we show that essentially the 
same bounds hold. Specifically, If m^ is the number of edges of 6k{S), then 

E[m|^]end|^/2±0(k^/n\ogn) , 

where d]^, defined above, is the average degree of the Q^-graph. in the Poisson model. We 
also give a concentration result that shows that the number of edges, m^t, is highly concen- 
trated around its expected value. In particular 

Fr{\m},-ndJ,/2\>k^/cn\ogn}<n~^'^''^ . 



1.2 Related Work 

As discussed in the introduction, a plethora of hterature exists on theta graphs and their 
applications, though most of this work focuses on worst-case analysis. One notable excep- 
tion is the work of Devroye et al. [13] who study the maximum degree of theta graphs and 
show that, if S is a set of n points independently and uniformly distributed in a certain 
unit square, then 0^(5) has maximum degree concentrated around ©((log«)/loglogM)Q 

In contrast, properties of other proximity graphs of random point sets have been 
studied extensively: 

• Devroye ['12] presents a general theorem for obtaining exact leading constants for 
the expected degree of a number of proximity graphs over point sets drawn from a 
large class of distributions. This theorem can be applied to Gabriel graphs, relative 
neighbourhood graphs, and nearest-neighbour graphs. This work [II 2[ Section 7] 
also mentions "directional nearest-neighbour graphs," now commonly known as Yao 
graphs [I15LI24J , and points out that the general theorem does not apply to these (nor 
does it apply to theta graphs — for the same reasons). 

• Penrose and Yukich [I21i 122] develop weak laws of large number and central limit 
theorems for several statistics of proximity graphs of random point sets under some 
assumptions about the locality of the graph and the statistic. Their results apply 
to statistics like total edge length and number of components of graphs such as k- 
nearest neighbour graphs, sphere of influence graphs, and Delaunay triangulations. 

• Bern et al. [6] study the maximum degree of Delaunay triangulations of random 
point sets. Devroye et al. study the maximum degree of Gabriel graphs [13] of ran- 
dom point sets. Arkin et al. [^J study the length of the longest edge in Delaunay 
triangulations of random point sets. 

• The issues of connectivity and giant components in the r-disk graph of random point 
sets — in which an edge uw is present if and only if the distance between u and w is at 
most r — is the subject of intensive research and there are at least two books devoted 
to the topic mi EHl. 

2 The Poisson Model 

In the Poisson model, the number of points in any region with area is A follows a Poisson 
distribution with parameter A. For definitions of point Poisson processes and distribu- 
tions see, for example, Daley and Vere-Jones [11 , Chapter 2]. For our purposes, the most 
important properties of the Poisson process are the following: 

1. The probability that a particular Lebesgue-measurable set, X, whose area (Lebesgue 
measure) is A is empty of points is exactly e^^. 



Devroye et al. actually consider the closely-related Yao graphs |15l I24| , but their proofs apply, almost 
without modification, to theta graphs. 




Figure 2: The edge uw is mutual if and only ii T'\T is empty of points. 

2. For disjoint regions Xi,...,Xy, the events "Xj is empty of points", for i e {l,...,v} are 
independent. 

Throughout this section, and in particular in Section [2^ we have made extensive 
use of Mathematica to do symbolic integration and manipulation of trigonometric func- 
tions. Mathematica notebooks containing the code for these calculations are available at 
the first author's webpage. 

2.1 Analysis of p/; for Even fc 

In this section we determine the value of pj. for even values of k. Somewhat surprisingly, 
the value of pi^ in this case does not depend on k. 

Lemma 1. For even integers fc > 4, p^ = ^ « 0.6045997883. 

Proof. Let u be an arbitrary vertex in a 0jt"gr^ph and let whe a vertex that u has chosen as 
a neighbour in one of its cones, C {u> is the "closest" vertex to u in C). Let T be the open 
isosceles triangle defined by C and a line through w that is orthogonal to the axis of C; 
refer to Figure [2] If the edge of T opposite u has length £, then w partitions this edge into 
two pieces of length r and £ -r. 

Let T' be the triangle obtained by reflecting T through the midpoint of the edge uw 
(so that w is a vertex of T'). Note that, because k is even, T' coincides with one of u's cones. 
In particular, if w e u + Ci, then T' is contained in w + Cj where j - {{i + k/2- 1) mod k) + l. 
Therefore, the edge uw is mutual if and only if T' \ T contains no points. The area of T' \ T 
is 

A{T'\T)^a{r^ + (£-rf) , 

where a = cos(0/2)/(4sin(0/2)). We now have enough information to compute the proba- 
bility that the edge uw is mutual conditional on £ and r: 



Vx{uw is mutual | (,, r\ 



:exp(-a(r2 + (^-r)2)) . 



Next, observe that the location of u> is uniformly distributed on the edge of T opposite u, 
so the value of r (conditioned on £) is uniform over [0,£], unconditioning r gives 



f{£) = Fr{uw is mutual | £} 

= (l/^)Pr{MU' is mutual I ^, r}dr 

Jo 



^0 

Jo 



(l/^)exp(-a(r2 + (^-r)2))dr 





= ^^ • exp(-a^V2) • erf(^V^) , 



where 



erf(A:) = -— e ^^ dz 
V^r Jo 

is the Gauss error function fTl Section 7.2]. 

Next, we remove the conditioning on £. The triangle T defines a region of area a£'^ 
that is empty of points. Therefore, by Property 1 of the Poisson process, the cumulative 
distribution function of £ is given by 

P(x) = Pr{£ < x} = 1 - exp{-ax^) , 

for X > 0. The probability density function of £ is therefore given by 

d 2 

p(x) = --P(x) = 2ax-exv)(-ax ) , 
ax 



for X > 0. Finally, we obtain p^ as 



Jo 
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Pk= \ p{£)-f{£)d£^— ^x 0.60 A5997 883 . D 



2.2 Analysis of p]^ for Odd k 

Next, we determine the values of pi^ for odd values of fc > 5. Although the strategy for 
doing this is the same as the even case, the odd case turns out to be considerably more 
complicated; the value of p;^ does, indeed depend on k. 



Lemma 2. For odd k >5, 



Pk = ^ 



where 



a 



' arctan(2(cos(f) + cos(f )f/(a^)) ] .„. 

+ arctan 4 2cos(^) + sin(^) \/{a^)\ ^^' 

7 = 4+11 cos I — I + cos I — I , 

/(27cos(f) + 17cos{f) + 3cos(f) + cos(^))csc{f) 



and 



VR¥J 



. /47r\ . /67r\ . /87r\ . llOnw 

-M8sm(-)+llsm(-) + sm(-) + sm( — )) . 



Proof. The proof proceeds in the same manner as the proof of Lemma [T] Let u be an 
arbitrary vertex in a 0/f-graph and let w be a vertex that u has chosen as a neighbour in 
one of its cones, C. Let T be the open isosceles triangle defined by C and a line through w 
that is orthogonal to the axis of C. See Figure [3] Assume that the side of T opposite u has 
length 2£. Using a suitable rotation, we may assume that the axis of C is horizontal and, 
by symmetry, we may assume that w is on or above the axis of C. 

Under the preceding assumptions, w is then uniformly distributed on a vertical 
segment of length € whose endpoints are on the axis of C and the upper boundary of C. 
Suppose that the distance from w to the axis of C is r. Then a straightforward, but tedious, 
calculation that mainly uses the law of sines shows that 

Vr[uw is mutual \€,r}- exp(-A - B) , 

where 

A — 

sin(27t/fc) 

and 

^= • rln^^l , •((^-r)cos(2n/fc) + ^)^ 
2sm(37r/fc)cos(7T//:) ^^ ' ^ ' ' 

This calculation is illustrated in Figure [Sland the accompanying worksheet shows the sim- 
plifications that lead to the expressions for A and B. Integrating over r gives us 



r2 



-i 



/(^) = Pr{MK; is mutual K}= (l/^)exp(-A-B)dr. (1) 

Jo 

Like the corresponding integral in the proof of Lemma [Tj (ITJ has a closed-form that in- 
cludes the Gauss error function. 

In order to remove the conditioning on (,, we need the probability density function 
for (,. Proceeding as before, we have the cumulative distribution function 

P(x) = PrK < x} = 1 -exp(A:Vtan(7r/fc)) , 



e-r 




Figure 3: The derivation of the Fr{uw is mutual | £, r} for odd k. 



and the probability density function 

Finally, we determine p^^ by integrating over €\ 

Pk- p{i)-f{i)d€ , 
Jo 

which (after introducing the variables a, ji, and y) yields the expression for pj^ given in the 
statement of the lemma. D 

3 Points in a Unit Square 

Next we argue that results similar to Lemmas [T] and [2| albeit with lower-order error terms, 
hold for the graph 0^(5), where S is a set of n points independently and uniformly dis- 
tributed in the square [0,y/n]^ of area n. Observe that, in this model, the probability that 
any particular region X c [0, V"]^ does not contain any points of S is exactly (1 -A/n)" - 
exp(-A) - 0{A/n), where A is the area of X. This is consistent with the Poisson model up 
to an additive error of 0{l/n). 

The primary work in this section involves finding quantities that look like those 
that appear in the previous section, but have an additive lower-order error term. To help 
manage these error terms, the notation x -y±a denotes that x is some value in the interval 
[y-a,y+a\. We will abuse this notation slightly by writing equations of the form x+a - y±b 
when [x -a,x + a]c[y -b,y + b]. Occasionally, we may also integrate expressions that use 

this notation. In this case, we use the inequality J (x±c)dx - \ x± c(b - a). 

We sometimes encounter expressions like A/(l - x), with < x < 1/2 which we 
bound by 

A/(l-x)^A + 0{Ax) . 

We also frequently encouter expressions like (1 -A/n)"~'^, where c is a constant and 
A < n/2. We will always bound these as follows: 

^ ' (I -AMY 

_ exp(-A)-0(A/n) 
~ {1-A/nY 
exp(-A)-0(A/n) 



_ exp(-A)-0(A/n) 

1 - 0(A/n) 
> exp(-A) - 0{A/n) 
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Figure 4: The support, [0, V"]^' is partitioned into a core, C, and a near-boundary C. A 
point u in the core almost surely has a neighbour in every cone, otherwise u is incident on 
a large triangle, T c [0, V"]^' that is empty of points. 



and, similarly. 



{\-Aln) 



< 



{\-Alnf 

{\-AlnY 

exp(-A) 

{l-AInf 

:exp(-A)(l + 0(A/n)) 



3.1 Expected Number of Edges 

In this section, we analyze the expected number of edges of Q]^[S). For each point u ^S 
and each i e {l,...,/c}, let e{u,i) be the edge (if any) that u chooses in its fth cone, u + Ci. 
That is, e{u, i) is the edge uw where w eu + Ci has the projection onto the axis of m + C, that 
is smallest among all points in S n u + Q. We define the height of the edge uw - e{u,i) as 
the distance between u and the of the orthogonal projection of w onto the axis of u + C;. 

For our analysis, we partition [0,^Jn]^ into a core, C = [2t, Vn-2t]^, where t = 
^Jcklogn, and a near-boundary, C = [0,^Jn]^ \ C (see Figure|4]|. The motivation for par- 
titioning into a core and near-boundary is that (1) there are not many points in the near- 
boundary and (2) points in the core behave almost exactly like points in the Poisson model. 
The following Lemma shows, for example, that points in the core nearly always have a 
neighbour in each of their cones. 

Lemma 3. For any u e S and any i e {l,...,k}, let £"„ ,■ denote the event "e{u,i) exists and has 
height at most t." Then Vr{£^^i | u 6 C} > 1 - m"^('^). 



Proof. Fix some location oi u e C and draw an open isosceles triangle, T, with apex u, 
contained in m + Cj, whose internal angle at u is = 2n/k and whose height is t (see 
Figured. Observe that, since u is in the core, T c [0, y/n]^. Furthermore, the area of T is 



A = t^tan{n/k) e @{t^/k) = e(clogn) 



For T to be empty of points in S, the n-\ points of S \ {«) must all fall outside of T. Thus, 
we have 



Pr{T nS=0|M6C} = (l - A/m)""1 

<exp(-A)(l + 0(A/n)) 

<2exp(-A) (for sufficiently large m) 

= 2exp(-0(clog«)) 



If £"„ f does not occur, this means that the event described above has occured. Therefore, 

Vx{£^^i I u 6 C} > 1 - «""<"', as required. D 



Our next lemma shows that edges generated by points in the core have essentially 
the same probability of being mutual as they do in the Poisson model. 



Lemma 4. For any u e S and any i 6 {!,..., fc), let Muj denote the event "£^,1 and e{u,i) is 
mutual." Then PrjA^,^ ,- | u 6 C} = pj. + 0{{\o^n)^/n). 



Proof. Throughout this proof, all probabilities we compute are conditional on u € C, even 
when this is not explicitly stated. The proof is basically a reproving of Lemmas [T] and 
[2] that takes care to deal with boundary effects. Here, we will prove the case for even k 
(Lemma [I]) only. The case for odd k (Lemma [2]) can be done the same way. 

We first compute the probability conditional on E^i and for fixed values of £ and 
r as described in the proof of Lemma [T] and illustrated in Figure [2] The notations £, r, T, 
and T' all have the same meaning as in the proof of LemmafXl The edge e{u, i) is mutual if 
and only if the remaining n-2 points in S \ {u,v;} — ^which are already conditioned on not 



10 



being in T — also fall outside of T'. The probability that this happens is 

,f!-2 



yr[Mu,r\£u,,^,r} 



l-area(TuT')/« 



1 -area(T)/« 



1 - area(T' \ T)/« - area(T)/« 



1- 



1 -area(T)/M 

area(T'\ T)/n^ 
1 -area(T)/n 



n-2 



n-2 



a(r^ + (£-rf)/n] 
1 - afVn 

a{r^ + {£-rf){l + 0{a£^/n)) 



n-2 



= exp(-a(r2 + (^ - rf){l + 0{afln))) + 0(a^V«) 
= exp(-a(r2 + (^ - r)^)) • exp(0(a2^VM)) + O(a^Vn) 
= exp(-a(r2 + (^ - r)2))(l + 0(a2^V«)) + 0{afln) 
= exp(-a(r2 + (^ - r)^)) + 0{{a^£^ + a^^y^) 
= exp(-a(r2 + {£- rf)) ± Oda^t"^ + at^)/n) . 



(2) 



(Step Q follows from the inequality 1 + (e - l)x > e^ for < x < 1.) We then remove the 
conditioning on r by integrating: 

f\£)^Vr{M,,\£u„e] 

= {l/£)(e^^{-a{r^ + {£ -rf)) + 0{{a^t^ + at^)/n))dr 

-J u 

= f[£) + 0[fya^t^ + at^)ln) , 

where j{£) is the same f{£) defined in the proof of Lemmafll 

To finish, we need the distribution function for £ conditional on <S„ ,-. The triangle 
T has area a£'^ so the probability that it is empty of points of S \ {u) is (1 - a£^h 
Therefore, for < x < t, we have the cumulative distribution function 



\n-\ 



P(x) = PrK<x|f„,,} 

_ Prjf",,^,- and £<x\ 

_ PrK<x} 

_ l-(l-axVn)""^ 
~ l-(l-at2/«)"-i 



(since 0<x<f, so^<x implies £"„;■) 
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From this we obtain the density function 

d 



p'(x).-P(x) 



2/^\n-l 



d 1 -(1 -ax In 
~dx I - (I - at^/n)"-'^ 
2ax{n- 1)(1 -axV«)"~'^ 

«(l-(l-at2/«)"~i) 
2ax{n - l)(exp(-ax^) ± 0(ax^/n)) 

2ax{n - l)(exp(-ax^) + 0{ax^/n)) 



n — n ^^y'-i 



2axM(exp(-ax ) + 0{ax /n)) 



Yl — Yl Li(c) 

.2\^rM — 2, 



(1-1/n) 



2ax(exp(-ax ) + 0(ax /«)) 



(1-1/n) 



1 _H-fi(c) 

= 2ax(exp(-ax2) + 0(axV«)) (1 + n""*'*)(l - 1/n) 

= 2ax(exp(-ax^) + 0((1 + ax^)/n)^ 

- 2axexp(-ax^) + 0((ax + a^x^)/n) 

= 2axexp(-ax^) + 0((at + a^t^)/n) 

^p{x) + 0(a^t^/n), 

where p(x) is the same p{x) that appears in the proof of Lemma [T] And now we have 
enough information to finish: 

Vr{Mu,i I u € C} = Vr{£,,^i \ueC}- Vr{Mu,i \ £u,i, ueC) 

+ Pr{4,,- I M e C} • Vr{Mu,i \ £u,u ueC) 
> Pr{£^^i \ueC]- Pr{Mu,i \ £u,u " e C} 
>il-n'^^'^)Vr{Mu,\£u,t,ueC} 

Jo 
,„-n(c)j r {^p^£) + o{a^t^/n))-{f{i) + 0{{a^t^ + at^)/n))d£ 

«-"W) r p(^)/(^) + O(p(0(a2t4 + af)/n + f{£)a^t'^/n + (a^t^ + ah^)/n^)d£ 
Jo 

'l -«-"<'') r p(^)/(^) + 0((a2t4 + at2y^^(^5j8^^3j5y„2j^^ 
7tV3 
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+ 0{{a^t^ + at^)/n) 



+ 0((logn)V«) 



12 




Figure 5: The quantities E^, Dq, and Mq: Eq counts the number of segments, Dq counts 
the number of arrowheads, and Mq counts the number of segments with two arrowheads. 



and 



Vr{Mu,r \ueC}< Pr{^,,, \ueC}- Vr{Mu,r \ £u,r] + (1 -Pr{^„,,}) 
< Pr{£„,, \ueC}- Vr{Mu, \ £u,] + """'"^ 

= ^±0((log«)V«) . 



D 



Lemma 5. Let S he a set of n points independently and uniformly distributed in [0, 1]^. Then 
the expected number of edges of 6^(5) is nd^l + 0(k^Jnk\ogn). 



Proof. In this proof, it will be helpful to distinguish between directed and undirected 
edges. Undirected edges are the edges 6]^{S) that we have been considering throughout. 
In contrast, uw - e{u,i) is a directed edge from u to w. If uw is mutual, then wu - e{w,i) is 
a distinct directed edge from w io u. In this way, if we let £ denote the number of undi- 
rected edges, D the number of directed edges, and M the number of mutual directed edges 
oi d^{S). Then we have 

E^D-M/2 . 



Let Eq, Dq, and Mq denote the same quantities but only counting those edges with at least 
one endpoint in the core. (See Figure [5]) 



13 



We begin with a lower bound: 

E[£]>E[£c] 

= E[Dc]-E[Mc]/2 



k 

= ^^Pr{" ^ C}{^A£u,i I u e C}-Vr{Mu,i \ u e C}/2) 

ueS i=l 
= k{^/n-Atf (?r{£,,j \ u e C}-Pr{Mu,i \ u e C}/2) 

> k(^/Jl-Atf(l - «-"("> - pk/2 - 0((log nfVn)) 

> {kn - 8t^fn)(l - n'^W - p,,/2 - 0((log «)t/«)) 
= kn{l-pk/2)-0{ty/n) 

- ndi(/2 - 0{^jnk\ogn) 

For the upper-bound we proceed as follows: 

E[£]<E[£c] + E[fc-|SnC|] 
^^[Ec] + k{St^fn-l6t^) 

< E[£c] + 4/c-y^McfcIogn 

< ndy2 + 0((log n) 2 ) + 4/c -y/nc/c log n 
= ndyt/2 + 0{k^Jnk\ogn) 

where the last step follows from a calculation similar to that done in the lower-bound. D 

3.2 Concentration 

Next, we show that the number of edges in this model is tightly concentrated about its ex- 
pected value. We begin with the following result, which follows immediately from Hoeffd- 
ing's Inequality [9\, and shows that the number of points in the core is highly concentrated 
around its expected value: 



Lemma 6. Pr|||S n C| - (V«-40^| >-^ck\o^n\en 



■n(c) 



To prove our concentration bound, we make use of a versatile concentration in- 
equality due to McDiarmid [1181 Lemma 1.2]: 

Lemma 7 (McDiarmid's Inequality). Let Xi,...,X„be independent all taking values in the set 
X and let f : X" i-^ M foe a function such that 

\f{xi,...,Xi,...,x„)-f{xi,...,x[,...,x„)\<k , 

for some k> and all Xi,...,x„,x^ e X and all f 6 {!,...,«}. Then, for all e > 0, 

Pr{|/(Xi,...,XJ-E[/(Xi,...,X„)]|>e}<2exp(-2eV(«/c2)) . 
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[o,v^^ 



Figure 6: The set Q ensures that every vertex in S has at least one point in each of its 
Oj^-cones. This figure shows a possible set Q for fc = 5. 

Lemma 8. Let S be a set of n points independently and uniformly distributed in [0, ^fnY and 
let nij^ denote the number of edges in d]^(S). Then Pr{|mj. -E[mj.]| > k^Jcn\ogn] < n~^^'^\ 

Proof. Let Q be a set of k points chosen such that, for any point u e [0, yfn]^, each of u's 
0jt-cones contains exactly one point in Q. (Q could be, for example, the vertices of a suffi- 
ciently large regular fc-gon. See Figurel6]) We begin by studying 6i^{S U Q). This graph is 
somewhat more nicely behaved since each vertex in S defines exactly k directed edges. 

For a set S = {ui,...,u„} C [0,y/n]^, let f(ui,...,u„) be the function that counts the 
number of edges in 6i^(S U Q). Observe that, for any u^, ...,u„, u'- e [0, 1]^, we have that 

|/(mi,...,m,-,...,u„)-/(mi,...,m-,...,m„)|</c . 

That is, moving Uj to location u'- can be though of as removing u,, resulting in the loss 
of at most k non-mutual edges emanating from u,, followed by adding u'- resulting in the 
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k 


Pfc 


dk 


dk/k 


4 


4 


0.60459978 


5.58160084 


1.39540021 


5.58190420 


5 


0.70168463 


6.49157685 


1.29831537 


6.49280635 


6 


0.60459978 


8.37240127 


1.39540021 


8.37587443 


7 


0.66778040 


9.32553719 


1.33221959 


9.32743266 


8 


0.60459978 


11.16320170 


1.39540021 


11.16338275 


9 


0.65740932 


12.08331611 


1.34259067 


12.08319579 


10 


0.60459978 


13.95400212 


1.39540021 


13.95657530 


11 


0.65259895 


14.82141149 


1.34740104 


14.82372518 


12 


0.60459978 


16.74480254 


1.39540021 


16.74624814 


13 


0.64993659 


17.55082423 


1.35006340 


17.55520514 


14 


0.60459978 


19.53560297 


1.39540021 


19.53693380 


15 


0.64830027 


20.27549589 


1.35169972 


20.27565544 


16 


0.60459978 


22.32640339 


1.39540021 


22.32956115 


17 


0.64722020 


22.99725644 


1.35277979 


22.99979580 


18 


0.60459978 


25.11720381 


1.39540021 


25.12062469 


19 


0.64646902 


25.71708858 


1.35353097 


25.71956002 


20 


0.60459978 


27.90800424 


1.39540021 


27.91343314 



Table 1: Numeric values of p^ and dj^. 

creation of at most k non-mutual edges emanating from u'-. Letting m^ denote the number 
of edges in dj^{S U Q), we immediately obtain, from McDiarmid's Inequality, 

Prj m^ -E[m^] > kyjcnlognl < 2exp(-2clogM) € «" ^^' . 



To extend this result to nij^, thenumber of edges of 0/;(S), we study Pr{|m^-mj.| > k^Jcn\ogn}. 
If |m^ - mj.| > k^cnlogn then the number of points of S in the near-boundary C exceeds 
-y/c«logn or 0j.(SUQ) contains edges that join points in Q to points in the core, C. Lemmap^ 
shows that the probability of the former event is at most m"^''^\ while Lemmaplshows that 
the probability of the latter event is at most n"^('^) D 

4 Discussion 

We have given exact closed-form expressions for the average degree of 0j."gr^phs for all 
k > A. It is easily shown that 0yt graphs with k = 1,2,3 are not spanners, so the cases 
fc > 4 are the most important. We have also shown that the number of edges in a dj^-graph 
of n points uniformly distributed in a square is highly concentrated around its expected 
value. These results can be used to inform practitioners about the optimal choice of k to 
use in a particular application. Table [T] gives the numerical values of dj^, pj^, and d^/k, for 
/:6{4,...,20}. 

Some of our results involved extensive symbolic manipulations that were done 
with the help of Mathematica. There is certainly the possibility of mistakes, either in the 
software or by its users. In order to help validate our analytical results, the final column 
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in TablefTlalso shows the results of the following experiment: 4,000 points were generated 
uniformly in the unit square [0, 1]^ and their 0j.-graph was computed. The average degree 
of points in the square [1/3,2/3]^ was then computed. The final column of Table [I] shows 
the average of these values taken over all 2,000 repetitions of this experiment. In all cases, 
it agrees with our theoretical results up to the first 3 digits. 

Yao graphs [1151 124] are closely related to theta-graphs. In the Yao graph, Y)((S), 
each vertex u is connected by an edge to the closest point in each of its theta cones, where 
closest is defined in terms of Euclidean distance. The same strategy we use to determine 
the average degree of 6^(5) can be applied to Yjt(S). Unfortunately, this fails to give closed 
form answers. In particular, when trying to repeat the proof of Lemma [Tj we obtain the 
formula 



poo r^t 

7[^ 0texp(-0tV2) 

Jo Jo 



p[^ \ 0texp(-0t72) exp(-dt^/2 + t^sin(y)sm{d-y)/sin{6))dydt , 
Jo Jo 

for which we are unable to obtain a closed form (here 6 - 2n/k). Nevertheless, one can 
evaluate this integral numerically to obtain estimates of the mutual edge probability and 
average degree in Yao graphs. 
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